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ACOUSTIC  RADIATION  FROM  FLUID  LOADED  RECTANGULAR  PLATES 


by 


Huw  G.  Davies 


Abstract 


The  acoustic  radiation  into  a  fluid  filled  infinite  half-space 
from  a  randomly  excited,  thin  rectangular  plate  inserted  in  an  infinite 
baffle  is  discussed.  The  analysis  is  based  on  the  in  vacuo  inodes  of 
the  plate.  The  modal  coupling  coefficients  are  evaluated  approximately 
at  both  low  and  high  (but  below  acoustic  critical)  frequencies.  An 
approximate  solution  of  the  resulting  infinite  set  of  linear  simul¬ 
taneous  equations  for  the  plate  modal  velocity  amplitudes  is  obtained 
in  terms  of  modal  admittances  of  the  plate-fluid  system.  These  admit¬ 
tances  describe  the  important  modal  coupling  due  to  both  fluid  inertia 
and  radiation  damping  effects.  The  effective  amount  of  coupling,  and 
hence  the  effective  radiation  damping  acting  on  a  mode,  depends  on 
the  relative  magnitudes  of  the  structural  damping,  i.e. ,  on  the  widths 
of  the  modal  resonance  peaks,  and  the  frequency  spacing  of  the  resonances. 
Expressions  are  obtained  for  the  spectral  density  of  the  radiated 
acoustic  power  for  the  particular  case  of  excitation  by  a  turbulent 
boundary  layer. 
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Introduction 


to  random  excitation 

Kuch  research  has  recently  been  done  on  the  response^of  struc¬ 
tures  vibrating  in  air.  In  general,  and  certainly  as  far  as  most 
practical  situations  are  concerned,  the  response  in  tnis  case  is 
effectively  the  in  vacuo  response,  no  consideration  of  the  interac¬ 
tion  of  the  structural  vibrations  with  the  associated  sound  'eld 
being  necessary.  The  acoustic  radiation,  if  required,  can  then  be 
estimated  from  the  already  determined  structural  response. 

Light  fluid  loading  effects  have  alsu  been  included  in  some  analyses. 
Much  of  this  work  is  based  on  the  statist  .cal  energy  method  of  Lyon 
and  Maidanik^.  Maidanik^  has  used  the  method  to  estimate  the 
radiation  from  finite  panels  vibrating  it  air.  More  recently, 

Leehey^  and  Davies  ^  have  discussed  its  application  to  turbulent 
boundary  layer  excited  panels.  The  present  analysis  is  to  a  certain 
extent  an  extension  of  some  of  this  work. 

Previous  work  concerning  water  loading  effects  has  concentrated 
on  spherical  sheila,  infinite  cylindrical  shells  and  infinite  thin 
plates  (see,  for  example,  Junger^,  which  contains  a  large  number  of 
additional  references, and  Maidanik^).  in  each  of  these  cases  the 
interaction  problem  is  simplified  because  the  in  vacuo  normal  modes 
of  the  structure  are  maintained  when  the  structure  is  submerged  in 
water,  as  the  acoustic  field  can  also  be  expanded  in  the  same  series 
of  character is tic  functions  or  modes.  For  the  infinite  plate,  the 
series  is,  of  course,  replaced  by  an  integral  over  a  continuous  spectrum 


of  wavenumbers. 
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In  situations,  such  as  that  to  be  considered  here,  whrre  the 
structural  response  is  characterised  by  a  discrete  wavenuab®’*  spec¬ 
trum  and  the  acoustic  field  by  a  continuous  wavenumber  spectrum,  the 
in  vacuo  normal  modes  are  not  retained.  However,  the  expansion  of 
fh<>  winery  r®«;>vn»®  «  structure  in  terms  of  its  in  vacuo  modes 

is  still  valid.  It  is  convenient  still  to  refer  to  these  functions 
as  modes  and  to  talk  of  the  resonance  frequencies  of  these  modes,  in 
which  case  although  we  do  not  refer  to  a  frequency  associated  with 
some  natural  mode  of  vibration,  we  still  imply  a  frequency  associated 
with  the  maximum  value  of  the  amplitude  response  of  a  characteristic 
function.  An  essential  feature  of  this  problem  now  becomes  the  coupling 
together  of  the  in  vacuo  modes  by  the  structure  fluid  interaction. 
Although  we  discuss  here  the  particular  case  of  radiation  into  a 
fluid  filled  semi- inf inite  space  from  a  rectangular  plate  in  an 
infinite  rigid  baffle,  the  arguments  given  concerning  the  amount  and 
the  effect  cf  the  modal  coupling  induced  can  obviously  be  applied  to 
other  geometries.  What  we  attempt  in  this  paper  is  hardly  a  complete 
solution  of  the  coupled  problem;  the  modal  interactions  are  far  too 
complicated;  but  rather,  after  a  considerable  but  necessary  series  of 
approximations,  we  determine  and  Interpret  the  most  important  fea¬ 
tures  of  the  9tructune-f luid  interaction  and  their  effect  on  the 
response  of  the  system. 

We  consider  a  simply  supported  thin  rectangular  plate  inserted 
in  an  infinite  rigid  baffle  and  fluid  loaded  on  one  side.  The  normal 
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vibration  velocity  field  of  the  plate  is  expanded  in  a  series  of 

the  in  vacuo  normal  modes  or  characteristic  functions  of  the  plate. 

This  approach  leads,  because  of  the  structure-fluid  interaction, 

to  an  infinite  set  of  simultaneous  linear  algebraic  equations  to  be 

solved  for  the  infinite  number  of  unknown  modal  response  amplitudes. 

These  equations  are  obtained  in  section  2,  below.  Furthermore,  many 

of  the  coefficients  in  these  equations;  those  .case  elated  with  the 

fluid  loading  terms;  are  defined  by  integrals  which  themselves  can 

only  be  evaluated  approximately  for  various  regimes  of  frequency. 

(2) 

Some  of  these  coefficients  have  been  evaluated  by  Maidanik  and 
(4) 

also  by  Davies  .  We  will  continue  to  refer  to  those  discussed  by 

(2) 

Maidanik  as  modal  radiation  coefficients  as  they  are  a  measure  of 
how  efficiently  a  particular  modal  shape  radiates  when  no  other  modes 
are  excited.  However,  we  also  require  the  modal  coupling  coef¬ 
ficients  connecting  the  vibration  of  one  plate  mode  with  that  of 
other  plate  modes  because  of  the  plate-fluid  interaction.  These 
additional  coefficients  are  obtained,  at  least  asymptotically  at 
low  and  high  frequencies,  and  discussed  in  section  3.  The  previously 
determined  modal  radiation  coefficients  can  obviously  be  obtained 
as  special  cases  of  the  modal  coupling  coef f icien«.u .  The  real  parts 
of  the  coefficients  are  associated  with  a  radiation  damping  effect 
on  the  plate  response:  the  imaginary  parts  lead  to  a  virtual  mass 
to  be  added  to  the  mass  of  the  plate,  hence  causing  a  decrease  in 
the  modal  resonance  frequencies. 
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In  sections  4  and  5  the  solution  of  the  resulting  infinite  set 
of  modal  equations  is  discussed:  for  low  frequencies  in  section  4, 
and  for  high  frequencies (but  below  the  acoutic  critical  frequency) 
in  section  5.  At  low  frequencies,  k^,  k  ^  <  n/2  (where  kQ  is 

the  acoustic  wavenumber  and  are  the  dimensions  of  the  plate), 

all  modes  have  similar  radiation  characteristics  and  the  modal  equa¬ 
tions  are  all  of  the  same  form.  An  approximate  solution  of  the  set 
of  equations  is  discussed  and  modal  admittances  for  the  plate- 
fluid  system  obtained  which  contain  the  important  coupling  effects. 
The  radiated  power  spectral  density  is  then  discussed.  A  simple 
expression  is  obtained  in  terms  of  the  modal  radiation  coefficients, 
the  modal  components  of  the  correlation  function  of  the  applied 
force  and  the  modal  admittance  functions.  It  is  assumed  that  the 
applied  force  is  such  that  it  causes  no  additional  modal  coupling. 
This  requires  that  the  typical  correlation  lengths  of  the  forcing 
field  be  much  less  than  the  panel  dimensions,  a  condition  that  is 
satisfied  in  many  practical  applications.  The  case  when  this  condi¬ 
tion  does  net  hold  is  discussed  briefly  in  an  appendix.  The  modal 
admittances  contain  the  virtual  mass  terms  and  additional  damping 
terms  due  to  the  fluid  loading.  The  inertia  coupling  terms  are 
small.  Because  of  the  coupling,  the  radiation  damping  term  is  found 
to  be  itself  a  summation  over  many  modes.  Now,  if  the  total  damping 
is  assumed  small,  the  power  radiated  in  a  narrow  band  of  frequencies 
is  mainly  from  the  modes  resonantly  excited  at  frequencies  within  the 
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band.  Furthermore,  we  need  only  consider  the  coupling  between 
resonant  modes  in  the  band.  It  Is  shown  that  the  magnitude  of  the 
radiation  damping  is  determined  by  the  amount  of  modal  coupling, 
that  is,  by  the  numbers  of  modes  that  interact.  This  in  turn 
depends  on  the  relative  magnitudes  of  the  structural  damping,  a 
measure  of  the  width  of  the  resonance  peaks,  and  the  frequency  spac¬ 
ing  between  resonance  peaks.  Under  the  assumption  of  light  damp¬ 
ing  there  can  be  no  likelihood  of  power  flow  between  modes,  that  is, 
no  coupling,  if  the  resonance  peaks  do  not  overlap.  This  dependence 
of  the  radiation  damping  on  the  structural  damping  is  discussed  in 
section  4.  Estimates  are  obtained  of  the  radiation  damping  under 
light  (hence  no  coupling)  and  heavy  structural  damping.  These 
are  then  used  in  the  expressions  for  the  spectral  density  of  radiated 
power  obtained  by  averaging  over  the  resonant  modes  in  narrow  fre¬ 
quency  bands . 

At  high  frequencies,  k  t  ,  k  l,  »  n,  the  modal  interactions 

o  I  o  J 

are  more  complicated,  and  an  analysis  of  the  modal  coupling  effects 

(2) 

correspondingly  more  difficult.  Maidanik  has  shown  that  at 
these  frequencies  the  modes  can  be  divided  into  three  groups  accord¬ 
ing  to  their  radiation  characteristics,  namely,  edge,  comer  and 
acoustically  fast  modes  (see  section  3,  below).  We  estimate,  in 
narrow  bands  of  frequency,  the  total  coupling  between  different 
types  of  modes.  Since  we  assume  small  damping,  it  is  sufficient  when 
considering  edge  and  corner  modes  in  the  estimates  to  include  onlv 
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resonant  modes.  We  restrict  the  analysis  to  frequencies  below  the 
acoustic  critical  frequency,  that  is,  the  resonant  modes  we  consider 
all  have  wavespeeds  on  the  plate  less  than  the  acoustic  wavespeed, 
and  thus  do  not  consider  the  case  of  resonantly  excited  acoustically 
fast  modes.  This  is  hardly  restriccive  in  practice;  the  acoustic 
critical  frequency  for  a  1/4"  steel  plate  in  water  is  about  400,000 
Herti.  However,  as  the  radiation  efficiency  of  the  acoustically  fast 
modes  is  high,  it  is  necessary  to  determine  whether  or  not  the  con.ri- 
bution  from  these  modes  to  the  radiated  power  at  any  frequency  is  of 
importance  even  when  the  modes  are  non-reeonantly  excited.  These  modes 
are  thus  included  in  the  analysis.  Our  estimates  oi  the  modal  coupling 
effects  in  narrow  frequency  bands  thus  include  the  interactions 
between  resonantly  excited  edge  and  comer  modes  and  non-resonantlv 
excited  acoustically  fast  modes.  We  find,  however,  that  in  most  cases 
the  acoustically  fast  modes  are  not  sufficiently  highly  excited  either 
by  their  being  coupled  to  resonant  edge  modes  or  by  the  acoustically 
fact  mode  component  of  the  external  forcing  field  to  give  any  considerable 
contribution  to  the  radiated  field.  Kxpressions  for  the  spectral 
density  of  radiated  power  are  again  obtained.  The  radiation  damping  of 
a  mode  is  again  found  in  the  form  of  a  summation  over  many  modes  because 
of  the  coupling  effect.  This  summation  is  valuated  tir  in  the  low 
frequency  case. 

In  section  4  the  directivity  of  the  radiated  field  is  brieflv 
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In  section  7  numerical  estimates  of  the  radiated  power  spectral 
density  are  made  assuming  the  plate  is  excited  by  a  turbulent  boundary 
layer.  Ccr^os^^  model  cf  the  correlation  function  of  the  wall  pres¬ 
sure  field  is  used.  The  estimates  are  compared  with  the  3pectral 
density  ootained  neglecting  fluid  loading.  A  wide  rar.ge  of  values  of 
the  structural  damping  is  used  to  demonstrate  the  dependence  of  the 
modal  radiation  damping  cu.  the  structural  damping. 

Finally  in  an  appendix,  we  discus3  the  case  whe~  the  modes  are 
also  coupled  together  by  the  external  field.  A  simple  expression  for 
the  radiated  power  spectral  uensity  can  still  be  obtained  at  low  fre¬ 
quencies  in  terms  of  the  modal  coupling  coefficients,  the  modal  admit- 

(g) 

tance  functions  already  discussed,  and  Powell's''  modal  joint  accep¬ 
tances.  It  would  seem  that  the  two  forms  of  coupling  act  independently, 
and  can  be  dealt  with  separately. 


The  Coupled  Equations  of  Motion 


A  simply  supported  thin  rectangular  panel  of  length  and  width 

is  inserted  in  a  flat  infinite  rigid  baffle.  The  acoustic  field 

* 

rsdioted  into  the  dense  fluid  in  the  semi-infinite  space  on  one  side 
of  the  panel  and  baffle  will  be  considered.  It  is  assumed  that  reither 
the  panel  vibration  nor  the  acoustic  field  affect  the  applied  external 

foi  zt. , 

Rr  'tangular  coordinates  (x^,  x^,  x^)  a  (x^)  are  chosen  as  usual, 
being  normal  to  the  panel  and  the  origin  being  at  one  comer  of  the 

* 

e .  g  water 


8 


panel.  The  fluid  fills  the  region  >  0  (see  Figure  1).  The  equation 
for  thin  panel  normal  displacement  y(x,t)  driven  by  a  pressure  field 
p(x,t)  is: 


v  Ml  a  “U 


\  (LSfc) 


T> 


U 


Z.l 


where  ?(x,  Xj,  t)  is  the  acoustic  pressure  generated  by  the  motion 
of  the  panel,  and  D  and  m^  are  the  flexural  rigidity  and  mass  per  unit 
area  of  the  panel,  respectively.  8  is  a  coefficient  introduced  to 
account  for  mechanical  damping  of  the  panel. 

Following  reference  A,  we  consider  the  frequency  Fourier  trans¬ 
form  of  equation  2.1  and  expand  the  panel  normax  velocity  displace¬ 
ment  v(x,uj)  in  terms  of  the  normalised  characteristic  functions 


M I  {  H  \ 

ij  T*v\  V  —  j 


ft; 


where  k  ■ 
m 


nm 


nTT 


,  k  m  j—  and  A  The  modal  equation  for  the 

^  «  I3  p  1  3  M 

frequency  Fourier  transform  of  panel  velocity  is  thus  obtained: 
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2  2  2 

where  k  «k  +  k  ,  n  is  a  modal  structural  loss  factor  and  p 

rnn  m  n  mn  -inn 

and  F  are  of  similar  form  and  defined  by 
nm 


T> 


P  can  be  expressed  in  terms  of  the  modal  velocity  amplitudes  v 
mn  7  r  nm 

using  the  wave  equation  for  P()c,  x^,  w)  in  the  acoustic  medium,  and  the 
boundary  condition  in  the  plane  of  the  plate,  namely 


V? 

b  K? 


O 


After  some  rearrangement,  we  can  obtain  P  (w)  in  the  form 

mn 


90 


T  ( »)  =  S  "R  .  r  (*)  V  ^  ^ 


l,r-i 


where  R  is  a  coupling  coefficient  connecting  the  (m,n)  mode  and 
mnqr 

the  (q,r)  mode  defined  by 


"R  (*',  --  — -  , 


oo 

rr 


swn(.b.)  s*  u) 


Z.5 


-  oo 
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Here,  Z(k^j)  Is  the  radiation  Impedance 


2( b,«)  *  p»c«  ( i  -  p.  } 


with  k  *|k|,  k  -  o)/c  and  S  (k)  is  a  shape  function  defined  by 
—  o  o  mn  — 


S  (w\ 


e-  ds, 


Z  4- 


Pr  .  .  ious  results  using  light  fluid  loading  approximations  have 
only  included  the  term  R  and  have  neglected  the  modal  coupling 
terms.  It  is  shown  below  that  many  of  the  modal  coupling  coefficients 
are  equal  in  magnitude  to  the  modal  radiation  coefficients. 

The  equations  of  motion  can  now  be  written  in  the  form 


|l)k  -  w <>»  -  4 f'V  Sdn *-i]  (►i)  ~  tK}  5^ ^  ^  ~  "*  ^  Kwi.  • 


In  the  following  sections  approximate  values  of  the  radiation 
coefficients  ate  obtained  for  both  high  and  low  frequencies,  and 

the  solution  of  the  system  of  coupled  equations  2.5  is  discussed. 
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3.  The  Radiation  Coefficients 

Values  of  the  modal  radiation  coefficients  have  been  obtained  in 

(2) 

references  2  and  k  for  various  frequency  regimes.  Maidanik's  classi¬ 
fication  of  the  modal  radiation  coefficients  is  shown  in  Figure  2,  The 
modal  wavenumbers  are  at  fixed  lattice  points  in  wavenumber  space. 

At  any  frequency  U)  we  define  the  acoustic  wavenumber  k^  *  ui/c^,  which, 

in  k-space,  divides  the  modes  into  those  with  k  <  k  ,  thus  having  a 
—  mn  o 

wavespeed  on  the  plate  greater  than  the  acoustic  wavespeed  and  hence 

called  acoustically  fast  modes,  and  those  with  k  >  k  called  acousti- 
J  mn  o 

cally  slow  modes.  The  acoustically  slow  modes  can  be  further  divided 
into  edge  and  comer  modes.  These  are  further  discussed  at  the  end  of 
this  section. 

The  types  of  approximations  used  in  references  2  and  k  are  also 

suitable  for  evaluating  many  'f  the  cross-coupling  coefficients.  However, 

(9) 

for  some  high  frequency  -esults  Lighthill's  methods  for  the  asymptotic 
evaluation  of  Fourier  transforms  of  generalised  functions  are  used. 

The  shape  fur  tlons  used  in  the  integral  for  Rnnqr  are: 


s.  (b)  = 


.  \  tk.jt 


The  Integrand  in  equation  2.3  thus  Includes  terms  of  the  form 


I  -v 


which  has  the  effective  value: 


2(1  -  (~l)ra  cos  k1£1) 
0 


if  m  and  q  are  either  both  odd 
or  both  even 

otherwise 


where  we  have  noted  that  the  rest  of  the  integrand  is  even  in  k^;  the 

sin  k,£,  term  in  the  integrand  thus  vanishing.  It  follows  that  R 

l  1  mnqr 

is  only  non-zero  when  both  m  and  q  and  n  and  r  have  the  same  parity , 
where  by  parity  we  mean  here  that  m  and  q(and  n  and  r)  must  be  either 
both  even  or  both  odd.  Each  mode  is  thus  coupled  to  at  most  only  one 
quarter  of  all  the  other  modes. 

The  radiation  coefficients  will  be  written: 


S 


mnqr 


mnqr 


+  i  T 

mnqr 
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CO 

nr 


ku.  P Qf„  k„*  k*  kc  k  ..  k 


r 


_('_' toS  Mi)(!~  (-i)rtwik^4)  dka  clk^ 

J)  w-  wXV-  K)K-  ki  Xk,'-  kr‘Xki-  V-kJ'* 


s  1 


The  real  part  of  the  integral,  8^  ,  involves  integrations  over 

values  |k|  <  k  .  We  note  that  the  only  singularity  of  the  integrand 
is  ihe  square  root  singularity.  Typical  values  of  the  function 


I  (kj) 
nv^  i! 


(tf-  K)(  k‘-  K‘) 


are  shown  In  Figure  3.  Th?  function  is  such  that 


CO 


1^0, )  A 


*  O 


W  f  \  • 


We  consider  the  term  S  first.  Various  approximations  are 

mnqr 

necessary  for  different  types  of  modes.  Evaluation  of  the  various 

coefficients  is  conveniently  divided  info  three  parts. 

a)  k  ,  k  »  k  (with  similar  results  applying  to  modes  with  k  ,  k 
n  r  o  m  q 

»  k  .)  k  i,  ,  k  i,  »  it  . 
o  o  1  o  3 

We  make  the  approximation  k*  -  k*  'v  -k*.  The  k^  integration 
can  then  easily  be  performed  leaving  the  approximate  result: 


&  paCe 


rvina  r 


k< 

kn  k  P 


I  (k.)  elk, 

»v\(y  V  / 


If  k  ,  k  <  k  we  can  write  the  Integral  as  /  -  /,  to  obtain  the 

in  c  o  °  ok 

o 

result: 


%  f0c0  k^  k<^  k0  f  J_ 

TT  A  v  k„  kr  ^  U-  $  k0 


where  6  is  a  Kronecker  delta.  The  last  term  here  will  be  neglected 
mq 

as  being  of  higher  order  in  pwers  of  (k^C) 

If  k  <  k  ,  k  >  k  ,  the  integral  is  written 
■  o  q  o 


1  1  ^ 


-  C£>ik(?,  ^  ! 


I 
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giving 


r 


&  Po  C  o  ^  ** 

TT  A  ^  kr  kq 


<v 


The 


case  ,  k  >  k^  is  treated  more  accurately  in  part  (c)  without 


the  restriction  k  JL  ,  k  9.  »  it 

o  1  o  “ 

b)  k  -  k  <  k  ,  R  JL  ,  k  1  »  7r 

m  .q  o  o  1  o  3 

Following  Maidanik^  ,  and  earlier  Kiaichnan^^  (his  equation 
5.5)  a  delta  function  approximation  can  be  used  in  the  form 


r*i  ry\ 


TT 
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leading  to  the  expression 
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For  k  ■  k  ,  a 
n  r 

tv.  the  acoustically 


similar  delta  function  approximation  leads  directly 

fast  mode  radiation  coefficient  S  -pc. 

nmmn  o  o 
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For  k  ,  k  >  k  we  obtain  the  edge  node  reault  of  equation  3.2, 
n  r  o 

For  the  coupling  coefficient  between  two  acoustically  fast  nodes 
or  an  edge  mode  and  an  acoutically  fast  mode  we  require  an  estimate 
of  the  Integral 
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COS  kjf*  )  3k; 


U/-io)  {w.*-  k/)"1 


for  k^  <  kQ.  As  this  integral  ia  relevant  only  for  edge  or  acoustically 

fast  modes,  it  is  sufficient  to  evaluate  it  for  large  values  of  k  l.,, 

o  3 

We  can  treat  the  terms  1  and  (-l)ncos  k^£^  separately  if  we  treat 
the  separate  integrals  as  Cauchy  principal  values.  The  first  integral 
is  zero  in  this  case  and  we  art  left  with  an  integral  which  we  can 


write,  following  Lighthill 


(9) 


as 
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where 
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and  H  Is  the  unit  step  function.  The  function  f  ha3  singularities 
at  zm  tk^/k^  and  r  *  1.  The  main  contributions  to  the  integral  at  large 

values  of  Can  be  estimated  from  the  Laurent  series  expansion  of 

f(r)  about  its  singularities.  We  find,  as  expected,  no  contribution 
from  the  singularities  at  r  -  t  k^/k^  .  and  from  the  square  root 
singularities  obtain  the  asvmptotic  result 


I  ^ 
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Equation  3. A  can  now  be  evaluated  for  all  cases.  The  result*  are 
summary  'd  below. 

c)  A  result  valid  for  all  comer  modes  at  all  frequencies  is 


obtained  from  equation  3.1  by  writing  it  in  the  form 
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We  note  that  because  of  the  parity  of  (m,n)  and  (q,r)  the  ’-1)  a  1 

(  i)n  terms  are  Interchangeable  with  ( —  1 ) Q  and  (-Pr  terras,  respectively 

S  is  thus  symmetric  in  ra  and  q  and  n  and  r. 

mnqr 

Values  of  the  real  narts  of  the  modal  'oupl  inc  coefficients  can 
be  snnmar ized  as  follows: 

If  i.ir. ,n)  is  an  N-tvoe  edge  mode  and... 
i)  iq.rl  is  an  edge  mode 
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ii)  (q,r)  is  a  corner  tn ode, 


2  po  Co  k 

TT^u  Urk, 


lii)  (q,r)  is  an  acoustically  fast  mode, 
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If  both  (m  n)  and  vq,r)  are  acoustically  fast  modes 


'V  < 


ii 

V  +  ('  OV  L  kr  kj,(  X  \  **<r«s  -  it's  SL. 

'  Vs;  WJ  '*h 


*  (-f  ^  kJJ  X.  \\os(Ur~)^ 

*At.  C  KeJ  *' 


2.  10 


-  - 1 


fM* 


r 


20 


Equation  (3.6)  gives  values  of  S  when  both  abides  are  corner 

em  qr 

modes  which  are  valid  at  el.’.  frequencies. 

We  note  that  the  edge, corner  and  acoustically  fast  modal  radiation 
coefficients  can  be  obtained  as  special  cases  of  the  above  results. 

Only  the  dominant  terms  in  the  express iona  for  different  modes 
have  been  retained.  At  this  point  we  make  no  comparisons  of  the 
relative  magnitude  of  the.  various  coupling  terms,  leavij  ■»  this  discus¬ 
sion  until  actions  4  and  5.  The  lelative  importance  of  the  coupling 
terms  will  obviously  be  influenced  by  the  numbers  of  each  of  the  various 
types  of  modes  which  are  coupled. 

In  a  similar  way  we  estimate  the  mass  loading  coupling  term 
Tmnqr*  ijl£e8ration  involved  is  over  all  wavenumbers  greater  than 

the  acoustic  wavenumber.  The  required  value  of  the  square  root  in 
equation  (3.1)  If  now  -i(k*  +  k*  ~  k*)  .  Thus ,  for  all  acoustically 

slow  modes  the  range  of  integration  includes  all  the  mod'l  wavenumbers. 
It  is  obvious  in  this  case  from  the  nature  of  the  integrand  that  the 
largest  coefficients  are  those  having  either,  or  both,  m*q  and  n*r. 

We  consider  first  the  coupling  coefficient  of  two  X~type  modes.  Hie 
range  of  integration  can  be.  divided  into  the  three  regions 
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For  two  X-type  edge  modes,  it  is  easily  seen  that  the  first  of  these 
terms  is  the  dominant  one.  W c  obtain  in  this  ease 
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We,  thus,  obtain  i  tegrals  of  a  form  similar  to  that  in  equation 
(3,5).  We  gain  treat  *-he  terms  1  and  (-1)°  cos  separately  as 

Gaud"  principal  values.  The  square  root  singularity  leads  to  the 
same  asymptotic  form  as  before.  However,  we  now  obtain  a  finite  contri¬ 
bution  from  the  term 

CO 

*• 

I 

-  d^5  ' 


TT' is  gives  the  dominant  term.  It:  is  sufficiently  accurately 


U" 


estimated  as 
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The  terms  arising  from  the  cos  integration  are  now  not  required. 

The  modal  mass  coupling  term  for  two  X-type  edge  modes  is  obtained 
in  the  form: 
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where  we  have  made  the  further  approximation 


MS 


We  will  be  concerned  in  the  fol'  wing  sections  with  the  coupling 
between  resonant  edge  modes  only,  that  is  with  modes  close  together  in 
wavenumber  space.  For  these  modes,  expression  (3.13)  is  approximately 
true. 

In  a  similar  way  we  may  obtain  the  mass  coupling  coefficients 
between  other  types  of  modes.  We  note  that  when  both  modes  are  corner 
modes  we  obtain  in  place  of  the  integral  in  (3,1.1,  ^proximate 

form 
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iii)  (q,r)  is  an  aeoutically  fast  mode 
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When  both  (ra,n)  and  (q.r)  are  comer  modes  we  obtain: 
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For  all  modes  we  also  have  the.  result: 


-  p*to 


rVul  kVtyA 


kiV>y\'^>  k0 


k  IVjy^  k.  (3 


The  mass  coupling  between  two  acoustically  fast  modes  is  negligibly 
small. 


2.  2° 


Again  as  in  the  expression  for  Sm  r  only  the  dominant  terms  in 

the  expressions  for  T  have  been  retained. 

ranqr 

Expressions  (3.6)  to  (3.10)  for  S  and  expressions  (3.16) 
r  mnqr 

to  (3.?0)  for  T  shov  a  marked  contrast  in  the  nature  of  the  coupling 
mnqr 

induced  by  the  radiation  terms  and  the  mass  loading  terms,  respectively. 
A  physical  basis  for  the  nature  of  the  radiation  coupling  coefficients 
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We  also  note  that  In  both  expressions  (3.11)  and  (3.14)  when  n*r  the 
delta  function  approximation  leads  directly  to  the  result 
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pac, 
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This  result  Is  true  for  acoustically  slow  modes.  Further  inspec¬ 
tion  3hows  that  is  negligibly  small  for  acoustically  fast  modes. 

Values  of  the  Imaginary  parts  of  the  modal  coupling  coefficients 
can  be  summarised  as  follows  (compare  equations  (3.7)  forward): 

If  (m,n)  is  an  X-type  edge  mode  and... 

1)  (q,r)  is  an  edge  mode 
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ii)  (q,r)  Is  a  comer  mode 
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can  be  given  in  terns  of  the  volume  velocity  cancellation  effects 

'2) 

described  by  Maidanik'  .  We  refer  to  Figure  4.  When  the  acoustic 
wavelength  A is  much  greater  than  the  modal  wavelength  A^  the  radia¬ 
tion  character  of  the.  mode  can  be  represented  by  a  series  of  pistons 

of  width  A  /2,  there  being  a  phase  change  of  180°  between  neighboring 
n 

pistons.  Afain  as  X  >  A  ,  the  combined  radiation  from  two  adjacent 
y  o  n  J 

half  pistons  is  of  dipole  character.  The  modal  radiation  character 
is  thus  described  by  a  series  of  dipoles,  leaving  at  each  end  mono¬ 
poles  of  width  A n/ 4 .  The  dominant  radiation  is  from  these  edge  monc- 
poles.  These  monopoles  are  themselves  coupled  if  A^  >  £^.  However, 

the  dipoles  further  combine  to  higher  order  .nultipoles  in  this  case 
and  the  dominant  radiation  is  again  from  the  edges.  If,  on  the  other 

hand,  A  <  A  ,  there  is  no  such  cancellation,  as  the  radiation  from 
o  n 

the  mode  snape  is  essentially  uncoupled  along  the  entire  length  of 

the  plate.  Ill  us  ,  tor  example ,  an  X-tyne  edge  mode  is  so  called 

(A  >  A  .  A  <  A  <  £_)  because  the  volume  velocity  cancellation 

In  the  direction  and  the  lack  of  cancellation  in  the  direction 

lead  to  dominant  radiation  from  strips  of  width  A  / 4  along  the  edges 

n 

x^=0  and  Xj=£^. 

We  mav  treat  the  modal  radiation  coupling  coefficient  >  as 

mnqr 

a  measure  of  the  contribution  from  the  (q,r)  modal  plate  velocity  to 
the  (m,n)  modal  component  of  the  acoustic  field.  W"  first  consider 
equation  (3.7),  where  both  modes  are  X-tvpe  edge  modes.  The  contribution 
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is  large  when  the  node  numbers  In  the  direction  are  the  same.  Along 
the  length  of  the  strip  radiators  we  require,  for  large  coupling, 
that  there  be  no  inter-modal  velocity  cancellations,  that  is,  the  strips 
must  vibrate  in  the  same  mode  shape  along  their  length.  Changes  in 
the  width  of  the  strip  are  of  negligible  importance.  Similar  arguments 
also  explain  the  presence  of  the  Kronecker  delta  functiors  in  equations 
(3.9)  and  (3.10)  .  However,  when  we  consider  the  radiation  coupling  of 
a  comer  mode  with  either  another  comer  mode  or  an  edge  mode  (equations 
(3.6)  and  (3.8),  respectively)  we  find,  of  course,  that  the  only  contri¬ 
bution  from  the  (q,r)  modal  plate  velocity  to  the  (m,n)  component  of 
the  acoustic  field  now  comes  from  the  uncancelled  roonopole  radiators 
in  each  comer.  The  only  criterion  for  there  to  be  a  coherent  contri¬ 
bution  is  that  these  comer  radiators  be  In  phase.  This  is  controlled 
solely  by  the  parity  of  (m,n)  a;.  ( q ,  t ) . 

The  imaginary  parts  of  the  modal  coupling  coefficients,  T  , 

mnqr 

lead  to  virtual  mass  terms  to  be  added  Co  the  mass  of  the  plate. 

Volume  velocity  cancellation  effects  are  now  of  no  consequence;  the 

inertia  terns  act  over  the  whole  area  of  the  plate.  This  is  demonstrated 

bv  die  fact  that  T  is  the  same  for  all  acoustically  slow  modes  and 

mnmn 

that  even  the  mass  coupling  terms  T  in  equations  (3. 16''  to  (3.1^1  are 

mnqr 

essentially  of  the  same  form,  irrespective  of  the  division  hv  radiation 
char. i  ter  1st  i' s  into  t  .g«.  and  corner  modes.  For  there  to  he  1  rge 
int  rtia  coupling  between  modes  we  require  that  two  mode  numbers  he 
the  same,  so  that  in  one  direction  the  modes  vibrate  in  the  same  shape. 

a  symbol  i  sod  bv  the  Kronecker  delta  functions  in  equations  3.16  to  3 . 1  *> . 

' 
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It  is  known  that  for  wave  motions  on  an  infinite  plate,  wave 
mo' ions  with  supersonic  phase  velocities  appear  as  a  damping  of  the 
plate  response  as  energy  can  be  transferred  from  the  plate  to  the 
fluid  becoming  acoustic  radiation,  whereas  wave  motions  with  subsonic 
phase  velocities  lead  to  a  virtual  mass  term  as  no  energy  is  trans¬ 
ferred.  This  observation  is  also  true  for  finite  plates,  the  addi¬ 
tional  constraint  of  the  edges  of  the  plate  being  sufficient  to  ■’btaiu 
acoustic  radiation  from  subsonic  wave  motions  on  the  plate.  The  volume- 

velocity  cancellations  that  lead  to  radiation  from  these  subsonic 

(2) 

wave  motions  have  been  described  above,  following  Maidanik  .  The 
essential  features  of  supersonic  and  subsonic  waves  on  the  plate  are 
still  preserved,  however,  as  shown  by  the  large  radiation  coefficient 
of  acousticallvfast  modes  (equation  3.5)  and  the  large  modal  inertia 
coefficient  of  acoustical;-/  slow  modes  (equation  3.20). 

The  forms  of  the  coupling  coefficients  suggest  that  the  set  of 
equations  (2.5)  can  be  conveniently  rewritten  in  the  form 
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The  solution  of  the  set  of  equations  (3.21)  is  discussed  in  the  fol¬ 
lowing  sections.  The  discussion  is  restricted  to  frequencies  below 
the  acoustic  critical  frequency. 


4 .  The  Low  Fre quency  Limit 


At  frequencies  such  that  k  SL  .  k  i.  <  tt  all  modes  are  of  comer 

o  1  o  3 

mode  radiation  character.  The  radiation  and  coupling  coefficients  are 
described  by  equations  (3.6)  and  (3.19).  Each  of  the  modal  equations 
(3. 21)  Is  thus  of  similar  fora.  They  can  be  written: 
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S  is  described  bv  equation  O.cO.  T‘  and  T  are  defined 
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by  the  expressions 


and 
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The  summations  in  equations  (4.1)  are  over  all  (q,r)  modes  with 

the  same  parity  as  (m,n) .  Thus,  there  are  four  separate  sets  of 

equations  (4.1)  to  be  solved  corresponding  to  the  four  combination- 

of  even  and  odd  values  of  m  and  n.  We  will,  of  course,  assume  that 

the  nomenclature  of  equation  (4.1">  covers  all  four  such  cases. 

The  set  of  equations  (4. 11  cannot  be  solved  exactlv.  We  ' .vk 

for  an  approximate  solution  of  these  equations  which  includes  the  ms-1 

features  of  the  coupline,.  The  solution  for  v  is  obviouslv  a  functio: 

mn 

of  all  the  modal  forcing  functions  n  .  However,  we  might  reasonably 

mn 

expect  that  the  solution  for  v  can  be  written  in  the  form  of  n  m-  '*1 

mn 

th 

.admittance  function  V  multiplving  the  (m,n!  modal  forclnc  function 

mn 

p  ,  with,  in  addition  higher  order  terms  of  sn.  Her  magnitude  invol- 
mn 

ving  the  coupling  coef f lc  ‘ents  and  ill  the  forcing  functions  p  .  We 
can  obtain  such  a  solution:  the  interpre tat  ion  of  the  approx  irat  ions 
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we  make  in  doing  so  is  more  fully  discussed  later  in  this  section. 
The  equations  (4.1)  are  first  rewritten  in  the  form 
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The  set  of  equations  (4.2)  can  be  rearranged  to  give  in 

terms  of  the  new  icrce  f  .  The  term 
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is  evaluated  bv  substituting  for  v  fro*,  the  equations  (4.2).  We 
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where  we  have  used  the  relationship 
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i  [It  may  be  helpful  to  note  here  that  the  sufficies  m,  q,  u  will  always 

[  refer  to  mode  numbers  in  the  x^  direction;  the  sufficie3  nr  r,  v  will 

always  refer  to  mode  numbers  in  the  direction,] 

If  we  now  substitute  from  equation  (4.3)  into  equation  (4.2) 
we  obtain  the  exact  form: 
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and  B  is  defined  by  the  expression 
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th»  fiad,  already,  an  ad.itt.nce  luncticn  containing  ,t  le„t 
°f  thC  Tmran  type  couPling  effect*. 

We  now  repeat  the  above  process.  The  tern 


r5  u~ 

is  evaluated  using  the  equations  (4.5).  When  we  substitute  for  this 
term  back  into  equation  (4.5),  noting  the  relation: 
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we  obtain,  the  exact  form: 
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and  B 


13 


is  defined  by  the  expression 
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Once  again  we  can  solve  the  equations 


(4.6)  as  is  symmetric 


and  satisfies  the  relation 
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We  finally  obtain  the  expression 
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The  equations  (4.7)  are  an  exact  relation  derived  from  equations 
(4.1)  However,  they  do  not  represent  a  solution  of  the  set  of  equa¬ 
tions  (4.1):  there  are  still  terms  in  v  on  the  right  hand  aide.  We 
could,  In  principle,  perhaps,  continue  "solving"  the  set  of  equations 
as  above.  However,  the  modal  velocity  term*  we  would  now  deal  with 
are  second  order,  that  is,  they  are  multiplied  by  products  of  the  form 

i 

S  T  .  We  would  thus  be  dealing  with  second  order  coupling 
ranqr  mnmv 

effects.  The  first  order  coupling  effects  can  be  extracted  directly 

from  equation  (4.7).  The  admittance  function  Y  we  assumed  could  be 
’  mn 

found  is  seen  to  be 
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The  remaining  terms  in  equation  (4.7)  include  the  excita¬ 

tion  of  the  model  velocity  v  by  all  the  modal  forcing  functions,  and 

mn 

also  tl.3  higher  order  coupling  effects. 

The  simplest  way  to  extract  fcb«  doml'-isat  terms  is  by  redoing 

the  above  analysis  and  discarding  the  higher  order  t  rms  as  the 

analysis  proceeds.  Thus,  in  equation  (4.5)  we  may  write  for  g _  the 

ion 

approximate  form 
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The  second  order  terns  we  neglect  here  are 
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From  the  definitions  of  the  coupling  terms  we  find  that 
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Later,  when  evaluating  the  magnitude  of  the  coupling  we  will  find 
it  sufficient  to  consider  only  the  coupling  between  resonant  modes.  Wa 
would  then  expect  that  for  resonant  mcdes  k^  and  in  equation  (420) 
are  of  similar  magnitude.  The  terms  we  neglect,  (4.9),  already  second 
order,  thus  also  tend  to  cancel  each  other. 

By  neglecting  other  similar  product  terms  we  can  finally  obtain 
an  approximate  form  of  equation  (4.7).  We  will  consider  this  to  be 
the  required  approximate  solution  of  the  set  cl  equations  (4.1),  This 
solution  is 
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and  %  is  defined  by  equation  (A. 8).  We  note  again  that  fhr  *  "*m  in  ex 
mn 

pression  3.20  is  not  included  in  the  definition  of  R 

mnmn 

Our  approximate  solution  is  in  terms  of  an  admittance  Y  -  -iwM« 

mn  '  mn 

which  itself  includes  both  modal  radiation  and  modal  inertia  coupling 

til 

effects.  We  have  also  retained  the  first  order  forces  on  the  (m,n) 
mode  of  the  other  modal  forcing  functions.  We  note  that  equation  (4.11) 
can  be  rewritten  in  the  form 
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However,  the  importance  of^  in  this  equation  is  obscured.  We  will 

mn 

find  when  considering  power  flow  into  and  out  of  the  plate  that  equa¬ 
tion  (A. 11)  is  a  more  useful  form.  The  summations  in  equation  (A. 11) 
are  concerned  with  the  forces  on  the  (m,n)  mode  due  to  the  modal  compo¬ 
nents  p  ,  and  the  forces  exerted  by  p  on  all  the  other  modes.  As 
fur  •» s  the  total  v.  iocity  field  of  the  plate  is  concerned  these  com¬ 
bined  effects  will  tend  to  be  self  cancelling.  We  then  find  that  the 
important  coupling  effects  are  contained  in'^L^.  This  becomes  evident 


i.en  we  consider  the  power  flow  in  the  system. 
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We  also  note  at  this  point  th**'  the  solution  (4.7)  can  be  obtained 
directly  by  evaluating 


Vr 

where  R  is  defined  by  equation  (4.12).  The  form  of  the  admittance 
mnqr 

function  does  not  become  as  easily  evident  in  this  case,  and  little  is 

saved  in  the  amount  of  algebra  required  as  we  must  still  expand 

1  3 

R  in  terms  of  S  ,  T  and  T  „  Having  obtained  the  quasi- 
mnqr  mnqr  mnmr  mnqn 

solution  (4.7),  however,  we  recognize  that  the  approximate  solution 
(4.11)  follows  by  assuming  that 
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The  error  terms  in  this  last  expression  are  just  those  terms, 
such  as  (4.10),  that  we  discarded  to  obtain  our  approximate  solution 
(4.11).  In  some  of  the  analysis  that  follows  we  will  make  direct 
use  of  the  approximate  relation  (4.14). 

Having  obtained  an  approximate  solution  for  the  modal  velocity 
field  of  the  plate,  we  now  use  this  solution  to  obtain  estimates  of  the 
total  power  input  to  the  system  and  the  radiated  acoustic  power.  For 
simplicity  in  obtaining  these  results  some  assumptions  will  be  made 
about  the  nature  of  the  forcing  field.  We  will  assume  that  it  is  both 
spatially  homogeneous  in  the  plane  of  the  plate  and  is  temporally  stationary, 
that  is,  its  correlation  function  is  a  function  of  both  spatial  and  tem- 
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poral  separation  only.  We  farther  assume  that  the  characteristic 
correlation  lengths,  including  the  corrected  length,  if  any,  in  which 
the  field  decays,  are  much  smaller  than  the  dimensions  of  the  plate. 

These  last  criteria  are  those  of  Dyer^\  These  assumptions  are  not 
unreasonable  and  correspond  well  to  typical  physical  situations.  They 
are  useful  in  that  they  lead  to  no  coupling  of  the  plate  modes  by  the 
forcing  field.  It  has  already  been  noted  that  we  assume  no  interaction 
between  the  external  force  and  either  the  plate  vibration  or  the  associated 
acoustic  field,  i.e. ,  the  external  force  is  specified  irrespective 
of  the  response  of  the  system. 

We  can  express  the  correlation  of  two  modal  forces  as  follows: 
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where  <X>  denotes  the  expected  value  of  X.  Because  of  our  assumption 
of  homogeneity  we  may  write 
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where  0(k,oj)  is  the  wavenumber- frequency  Fourier  transform  of  the  cor¬ 
relation  function  of  the  force  p(x,t).  Dyer's  criteria  imply 
that  <p(k,ui)  is  essentially  constant  in  k.  The  orthogonality  of  the 
characteristic  functions  then  leads  to  the  simple  result 

-  tv-  tr  > 

where 
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f  (b,«)  |  | 

is  a  modal  correlation  function. 

*  ( 8  ) 

The  correlation  functions  <p  pT  >  are  essentially  Powell's 

mn  qt 

joint  acceptances.  When  Dyer's  criteria  are  not  satisfied  we  must 
include  all  these  functions  in  the  analysis.  However,  as  these 
criteria  are  satisfied  in  a  large  range  f  physical  situations  it  is 
assumed  ir  the  analysis  that  follows  that  the  joint  acceptances  always 
satisfy  the  simple  relation  (4.15).  We  may,  indeed,  reasonably  conjec¬ 
ture  that  the  two  modal  coupling  mechanisms,  the  external  field  and 
the  acoustic  field,  act  independently,  at  least  to  first  order.  Thus, 
as  far  as  a  study  of  fluid  loading  is  concerned  there  is  no  loss  of 
generality  in  assuming  that  Over's  criteria  are  met.  We  expect  the 
extension  of  the  analysis  to  include  the  additional  coupling  effects 


AO 


to  be  straightforward,  although,  possibly,  laborious.  An  approximate 
solution  is  suggested  in  an  appendix  for  cases  when  Dyer's  criteria 
are  not  satisfied. 

We  now  estimate  the  energy  flow  into  and  out  of  the  plate.  We 
consider  first  the  radiated  acoustic  power.  The  total  power  flow  into 
the  fluid  from  tne  plate  is  defined  to  be  an  integral  over  the  area 
of  the  pla^e,  namely 
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For  the  present  we  will  not  restrict  attention  to  the  real  part  of  this 
expression  but  rather  consider  the  total  complex  power.  We  may  refer 
to  the  real  and  imaginary  components  as  the  resistive  power  and  reactive 
power,  respectively.  The  resistive  power  describes  the  drain  of 
energy  by  radiation  to  infinity;  the  reactive  power  is  a  measure  of 
the  energy  that  is  stored  in  the  fluid.  To  obtain  the  spectral  density 
of  the  radiated  power, H  (u) ,  we  require  the  expected  value  of  P(x,  0,  t) 
times  v(x,  t  +  l)  for  all  time  delays  T.  We  again  consider  both  the 
real  and  imaginary  parts  of  H(ai). 
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By  using  the  modal  expressions  and  the  orthogonal i ty  ol  the 
characteristic  functions  we  obtain  the  result 
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where  the  summation  is  over  all  modes 
is  given  by  the  expression 
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P  in  our  present  notation 
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where  R  is  defined  bv  equation  (4.12)  and  we  have  ■’ociuded  the 
mnqr 

term  T  separately.  If  we  substitute  for  v  using  equation  (4.11) 
mnmn  qr 

we  obtain  directly 


S' 


\  r 
.i. 


i  r 
\ 


U>  now  use  the  fact  that  our  solution  is  based  on  the  approximate 
relatin'  (4.14).  To  this  degree  of  approximation  the  double  sin  at  ion 
above  can  be  reelected.  A  similar  double  summation  is  neglected  when  we 
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substitute  for  P  and  v  in  equation  (4.16).  When  Dver's  crit  i 
mn  mn 

are  assumed  to  hold,  the  expression  for  the  (complex)  radiated  power 
reduces  to  the  simple  form 
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The  real  part  of  this  expression,  the  spectral  density  o’;  the  power 
radiated  to  infinity,  say  n(cn)  ,  is  seen  to  he 


v..r  4 


where,  following  I.eehev  ‘  ,  we  have  written  ?  m  r  c  c  ,0  being 

mnmn  o  o  mn  mn 

the  modal  radiation  efficiency  (n  -  1  for  acoustically  fast  modes'1. 

.nr 

Kquat  ion  (4,18)  is  similar  in  form;  to  previously  ov,'Mshed  results  for 
the  radiated  power  obtained  either  by  neglecting  fluid  load  *nj>  effects 
or  bv  assuming  light  fluid  loading.  The  importance  of  the  admittance 
funct  ion  i;  now  clearlv  see-..  The  fluid  loading  effects,  the  modal  c  cur¬ 
ling  indv.ved  both  bv  radiation  loading  and  inertia  loading,  are  all 
included  in  the  admittance  functions  Y 

mn 

For  a  specified  external  force  acting  on  the  plate  the  spectral 


dens  i  *  v  o!  a  roust  ic  power  radiated  bv  the  plate  can  he  ciw.ruted  directly 


from  equation  (4.18).  When  t'.f  plate  modal  density  in  frequency  space 
is  high,  more  useful  express  io  is  can  be  obtained  by  averaging  IT(u))  over 
narrow  bands  of  frequency  and  considering  only  the  power  radiated  by 
modes  with  resonanace  frequencies  within  tne  frequency  band.  It  is 
assumed  th<_t  even  with  the  sided  damping  due  to  fluid  loading  the  resonance 
peaks  are  sufficiently  marked  that  the  radiated  power  is  indeed  mainly 
from  the  resonant  modes.  It  is  again  emphasized  that  the  resonance 
frequencies  referred  to  are  merely  the  frequencies  at  which  the  real 
components  of  the  admittance  functions  Y  vanish.  The  inertia  terms, 
of  course,  cause  a  decrease  this  frequency  from  the  in  vacuo 
resonance  frequency.  By  averaging  over  a  narrow  band  of  frequencies,  .(tu , 
equation  (4.18)  for  the  radiated  power  spectral  density  can  be  written 


where  the  summation  is  over  those  modes  with  resonance  frequencies  in  the 
Aoj  band. 
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The  integration  over  jY^|  will  now  be  carried  out.  To  emphasize 

the  coupling  effect  we  write  Y  in  the  form 

mn 
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where  T  is  replaced  by  p  c  T  .  The  terms  p  c  T 

qrqr  o  o  qr  o  o  qr 


are,  in  general ,  con 


siderably  smaller  than  tne  inertia  term  p  c  t —  included  in  B  .  The 

o  o  k  an 

mn 

effect  of  inertial  coupling  will  thus,  presumably,  be  small.  When 

averaging,  as  in  equation  (4.19),  we  are  interested  in  those  modes  with 

resonance  frequencies  in  &ij.  Conversely,  when  estimating  the  coupling, 

as  in  equation  (4.20)  we  are  interested  in  frequencies  near  the  resonance 

frequency  of  Y  ,  and  hence,  as T  is  small,  near  the  resonance 

freouency  of  B  .  The  added  damping  is  seen  to  be  due  too  ,  as  in 
mn  mn 

light  fluid  loading  approximations,  plus  other  modal  radiation  coef- 

B 

ficientscr  modified  by  the  terms  To  evaluate  the  expression 

qr  B 

qr 

(4.20),  the  frequency  is  written  in  terms  of  the  resonance  frequency 

uj  of  B  , 
mn  mn 


Pi  =  ry. 
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where  e  corresponds  to  the  small  variations  of  frequency  within  the 

do  band.  Hie  resonance  frequency  of  B  is  related  to  that  of  » 

qr  mn 

by  the  expression 
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where  X  is  a  measure  of  the  separation  of  the  resonance  peaks.  At 
frequencies  within  Au  we  can  now  write 
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and,  similarly, 
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The  summation  in  equation  (4.20)  can  now  be  expressed  in  the  form 
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It  is  clearly  ssen  that  the  important  contributions  to  the  summation 

at  the  resonance  frequency  of  B  come  from  those  modes  with  resonance  ' 

ran 

frequencies  close  to  Thus,  the  amount  of  coupling  is  controlled 

by  the  spacing  of  the  resonance  frequencies,  there  being  a  large  coup¬ 
ling  effect  into  those  modes  with  A™  <<  n  ,  that  is,  when  the  separa- 
°  qr  mn  ’  r 

tion  is  less  than  the  effective  widths  of  the  resonance  peaks.  On  the 
other  hand,  since  we  assume  that  the  modes  are  only  appreciably  excited 
at  their  resonance  frequencies,  we  neglect  the  coupling  between  modes 
whose  resonance  peaks  do  not  overlap. 

Two  cases  must  be  distinguished  in  evaluating  expression  (4.25). 


When  the  structural  damping  is  small,  that  is  when  A1"11  »  n  for  all 

qr  mn 
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inodes,  then  there  la  no  effective  modal  coupling,  and  expression  (4.25) 
reduces  to  a  single  term.  We  obtain  directly  in  this  case  that  at 
frequencies  near  resonance: 
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where  we  have  dropped  the  tern  p  c  X  as  it  la  much  smaller  jhan  the 
,  r  o  o  ran 

k 

term  p  c  •; —  .to  is ,  of  course,  already  defined  with  the  inertia 
o  o  k  ran 

ran 

k 

term  p  c  —  included.  Expression  (4.26)  is  not  strictly  complete  as 
o  o  k 
mn 

we  have  not  included  the  possibility  o:  two  modes  having  the  same 
resonance  frequency,  in  particular  as  in  the  degenerate  case  of  a 
square  plate.  However,  at  very  low  values  of  the  structural  damping 
we  see  that  we  may  reasonably  neglect  modal  coupling.  The  particular 
case  of  a  square  plate  can  be  included  in  this  analysis  simply  by  doub¬ 
ling  the  value  of  the  radiation  damping. 

The  second  case  concerns  plates  with  large  structural  damping  and 
high  modal  densities.  Considerable  modal  interaction  now  occurs.  To 
evaluate  the  coupling  effects  we  may  suppose  the  modes  evenly  spaced 
in  frequency,  in  which  case  the  separation  between  adjacent  modes  is 
obtained  from  equation  (4.22)  as 
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where  y  is  the  frequency  modal  density. 

A  »  n  ,  that  is ,  we  have  the  condition 
irm 

V  »  l 


We  are  now  interested  in  the  case 

Z ? 


The  summation  (4.25)  can  be  treated  approximately  as  an  integral 
in  A  in  the  form 
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The  limits  have  been  taken  as  ±  00  to  include  all  modes,  a  factor 

1/4  being  introduced  to  account  for  the  fact  that  there  is  coupling 

only  between  modes  of  the  same  parity.  The  imaginary  part  of  the  integral 

involves  the  factor  (  »  -  n  )  and  has  been  neglected.  We  thus  obtain 

inn  qr 

for  the  admittance  function  near  the  resonance  frequency 


Again  we  hove  drooped  Che  additional  inertia  terns,  r  .  Had  ve 
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Included  these  terms,  the  effective  mass  of  the  plate  would  be  described 

by 
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Typically  we  expect  that  the  last  term  here  is  small,  but,  if  necessary, 

it  can  be  included  in  the  estimate  of  the  resonance  frequency  tc  .  Equa- 

rnn 

tion  (4.28)  still  correctly  describee  the  nature  t  the  function  Y  near 

ran 

its  resonance  frequency  although  this  resonance  frequency  is  now  slightly 

changed.  The  important  fluid  loading  inertia  term,  fhe  modal  "self" 
pc  k 

inertia  — —  ,  has  already  been  included, 

mn 

The  dominant  fluid  loading  effects  are  the  change  of  resonance 

frequency  by  the  modal  "self"  inertia  and  the,  possibly,  large  increase 

in  the  total  damping.  Some  typical  magnitudes  of  these  effects  e"e 

discussed  In  section  7.  The  inertia  modal  coupling  terms,  although 

adding  considerable  complexity  to  the  analysis,  are  found  to  have  little 

effect  on  the  admittance  of  the  plate-fluid  system.  In  this  context 

1  3 

it  Is  worth  noting  that  if  we  neglect  the  terms  T  and  T 

°  °  mnmr  mnqr 

throughout  the  analysis,  equations  (  4.8),  (4.11)  and  (4,12)  then 

together  represent  an  exact  solution  of  the  set  of  modal  equations  (4.1). 

2 

The  integration  over  |Y  |  (now  as  a  function  of  e)  can  now  be 


performed  as  in  previous  analyses  (e.g.  references  3  and  4)  by  assuming 
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that  the  total  damping  is  sufficiently  small  that  the  effective  reso¬ 
nant  bandwidth  is  much  less  than  the  bandwidth  Aw.  We  obtain  two  estimates 
of  the  radiated  power  spectral  density  by  using  expressions  (4.26)  and 
(4.28),  respectively: 
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The  summations  are  over  all  the  resonant  modes  in  the  Aw  frequency  band. 

In  a  similar  way  we  can  obtain  an  estimate  of  the  total  input 
power  to  the  plate-fluid  system.  Thus,  corresponding  to  equation  (4.16): 
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The  required  power  is  the  real  part.  This  expression  can  rlso 
be  est  mat'  1  by  averaging  over  the  frequency  band.  The  Integral 
over  the  second  term  is  ter-  .  We  obtain  an  expression  for  the  spectral 
density  of  input  power  which  is  independent  of  either  the  structural  or 
the  radiation  damping:  the  only  fluid  loading  effect  is  the  change  of 
nodal  resonance  frequency  by  the  modal  self  inertia  term.  Thus, 
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It  is  shown  in  reference  4  that  the  modal  correlation  functions 

$  (uj)  re  the  same  for  all  resonant  modes  in  a  narrow  frequency  b_nd 
mn 

if  the  correlation  lengths  of  the  forcing  field  are  considerably  less 

than  the  corresponding  modal  wavelengths.  Fox  excitation  by  a  turbulent 

boundary  layer  tnis  condition  is  satisfied  at  frequencies  well  above 

the  hydrodynamic  critical  frequency,  that  is,  the  resonance  frequency 

of  tne  plate  modes  with  wavespeed  on  the  plat^  equal  to  the  convection 

velocity  of  the  forcing  field.  In  this  case,  the  functions^  can  be 

mn 

taken  outside  the  summations  in  equations  (4.24)  and  (4.25).  The 
radiated  power  spectral  density  is  then  linearly  related  to  the  input 
power  soectra"!  density. 

A  comparison  of  some  typical  numerical  values  of  the  radiated  power 
spectral  density  showing  the  magnitude  of  the  fluid  loading  terms  is 
made  in  section  7. 
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5 .  The  High  Frequence  Limit 

At  high  frequencies,  kQ£.p  k^f.^  >>  ti  ,  the  radiation  and  coupling 
characteristics  of  tue  inodes  are  not  the  same  for  all  modes  as  ui  the 
low  frequency  case.  The  coupled  mo^al  equations  are  thus  considerably 
more  complex  and  cc  “respond ingly  more  difficult  to  solve.  In  this  sec¬ 
tion  an  approximate  solution  of  ne  modal  equations  we  have  obtained  >' 
discussed.  We  start  by  est  iat  ng  the  magnitude  of  the  force  axer  ad 
on  a  mode  due  to  its  being  couf ie  to  other  types  of  modes.  Thus  t  r 
example,  we  make  an  estimate  of  t  ,e  total  affect  on  the  response  of 

an  edge  mode  of  the  coupling  into  all  the  corner  modes  ,  based  on  a 

typical  comer  mode  resonant  amplitude  and  the  total  number  ?f  resonant 
coupled  comer  modes.  This  approach  leads  to  a  first  approximation  in 
which  the  comer  and  edge  mode  responses  can  be  solved  for  s.  aratelv, 
including  cnly  the  coupling  into  other  comer  and  edge  roues,  respectively. 
We  include  the  nonresonant  acoustically  fast  modes  in  the  analysis  and 
calculate  their  response  due  to  coupling  to  resonant  edge  mode. .  This 
step  seemed  necessary  as  the  associated  radiation  from  these  highly 
efficiently  rad iat ig  modes  could  be  of  importance.  An  estimate  of  the 

total  radiated  power  shows,  however,  that  the  effect  of  these  modes  Is 

negligible.  From  some  numerical  estimates  o'  the  radiated  powei  made 
in  section  7  it  would  also  seem  that  the  modal  component  of  the  force 
exciting  the  nonresonant  acoustically  fast  modes  also  leads  to  a  neg¬ 
ligibly  smell  contribution  to  the  power. 
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The  comer  mode  contribution  to  the  radiated  power  is  included 
in  this  high  frequency  analysis.  It  is  known  that  under  light  fluid 
loading  conditions  at  frequencies  where  edge  mode  radiation  occurs, 
this  radiation  dominates  that  from  comer  modes.  Under  dense  fluid 
loading  conditions,  however,  the  higher  radiating  efficiency  of  the 
edge  modes  causes  correspondingly  higher  damping  of  the  modes,  the 
effect  being  such  that,  in  some  cases,  the  edge  mode  radiation  domimance 
does  not  occur  except  at  frequencies  where  very  many  edge  modes  are 
excited.  Structural  damping  is  important  in  determining  the  relative 
magnitudes  of  the  edge  and  comer  mode  radiation.  This  is  discussed 
in  section  7. 

We  first  consider  the  effect  on  the  edge  mode  response  of  the 
coupling  into  other  tvpes  of  modes.  The  complete  modal  equation  for 
an  X-type  edge  mode  response  will  be  writter 


t- 


fV* 


)( 

Here,  v  is  the  modal  velocitv  amplitude  of  an  X-tvpe  edge  mode, 
mn 

X  F 

R  is  the  coupling  coefficient  between  the  (m,n)  X-tvne  edge  mode 

mnqr 

and  the  (q,rl  acoustically  fast  rode:  similar  interpretations  follow 
obviously  for  the  other  notations,  C  Indicating  a  comer  mode.  The 
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summations  are  over  all  the  coupled  modes  taking  the  parity  of  the  mode 
numbers  into  account.  Similar  equations  exist  for  Y-type  edge  modes, 
comer  modes  and  acoustically  fast  modes. 

Estimates  of  the  magnitudes  of  the  summations  inequation  (5.11  are 
made  using  the  values  of  the  coupling  coefficients  obtained  in  section  3. 
Thus,  using  equations  (3.7)  and  (3.16), 
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using  equations  (3,8)  and  (1.17), 
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and  using  equations  (3.°)  and  (1,18) 


; 

* 


l’  ^ .  .  v 


r- 


54 


Now,  if  there  is  a  considerable  amount  of  modal  interaction,  the 

relative  importance  of  the  terms  in  expressions  (5.2),  (5.3)  and  (5.4) 

can  be  estimated  by  assuming  that  all  the  acoustically  slow  modes  in 

a  narrow  wavenumber  band  Ak  are  coupled  together.  We  also  include  all 

* 

the  nonresonant  acoustically  fast  modes.  The  case  of  ?ow  structural 
damping  with  essentially  no  modal  interaction  can  then  be  obtained  a a 
a  particular  case.  We  assume  when  estimating  the  relative  importance 
of  the  various  terms  that  all  the  resonant  acoustically  slow  moaes  and 
all  the  non-resonant  acoustically  fast  modes  have  the  same  velocity 
amplitude,  respectively.  The  radiation  coupling  and  the  inertia  coupling 
are  estimated  separately.  The  following  ratios  are  obtained  from 
equations  (5.25 ,  (5.3)  and  (5.4): 


i)  radiation  coupling  force  on  X-tvpe  edge  mode 
•*  C  ^ 
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inert .a  coupling  force  on  X-tvpe  edge  mode 


We  cannot  yet  estimate  the  effect  of  the  acoustically  fast  modes 


on  the  edge  ...ode  response  without  first  considering  the  acoustically 
fast  mode  equation.  This  is: 
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The  acoustically  fast  modes  we  treat  are  at  frequencies  well 

above  their  resonance  frequencies:  they  are  thus  mass  controlled.  This 

2 

is  reflected  in  the  term  (u  m  +  iup  c  )  which  Is  merely  the  mass  con- 

poo 

trolled  admittance  plus  the  modal  radiation  term  cl-p  c  .  'r*  ere  is  no 

oo 

virtual  mass  term  for  these  modes  The  summations  in  equation  (5.7) 
can  be  estimated  as  before.  However,  the  dominant  term  in  this  equa¬ 


tion  is  the  modal  admittance:  a  comparison  of  the  terms  shows  that  v 

-1  X 

is  et  least  a  factor  (k  £)  times  v  .  That  the  nonr»3npar't  amplitudes 

O  HiH. 

F  % 

v  are  very  much  smaller  than  the  resonant  amplitudes  of  v  is  hardlv 
qr  mn 

a  surprising  fact.  But  it  makes  possible  a  number  of  approximations. 

F 


In  equation  (5.7),  the  summation  E  is  dropped.  Only  the  inertia  coup- 


X  Y 

ling  terms  are  retained  in  E  and  Er  the  radiation  coupling  terms  are 

negligible  compared  with  the  radiation  term  up  c  .  The  edge  mode 

o  o 

inertia  coupling  is  retained  ss  being  of  possible  importance,  particu¬ 
larly  if  the  modal  force  p^  is  small.  Returning  to  equations  (5.5) 

qr 


and  (5.b),  we  now  find  that  where  there  ia  considerable  modal  Interaction, 

so  that  Ak  is  reasonably  large,  the  edge  mode  -^nations  dominate, 

although  not  to  any  great  extent.  For  example,  if  ■  k^/2  find 

k  -  10  k  ,  the  ratios  of  edge  to  corner  mode  summations  in  equations 

(5.5)  and  (5.6)  are,  respectively  1  :  1/10  and  1  :  k^/1200  Ak.  It 

is  thus  reasonable  to  Include  only  the  edge  mode  summation  in  the  edge 

mode  equation  (5.1).  (As  we  are  considering  frequencies  such  that 

k  £,  ,  k  £_  »  7T  the  acoustically  fast  mode  summations  are  obviously 

o  1  o  3  3  3 

negligible.)  In  a  similar  way  we  find  that  the  dominant  coupling  in 

the  modal  equation  for  the  corner  mode  response  is  into  other  comer 

modes.  We  have  thus  reduce’  the  problem  at  high  frequencies  to  solving 
for  the  edge  and  comer  mode  responses  separately,  neglecting  the  coup¬ 
ling  between  edge  and  comer  modes.  Having  found  the  edge  mode  response 
we  can  then  estimate  the  nonresonant  acoustically  fast  mode  response. 

The  acoustically  fast  modes  are  retained  in  the  analysis  until  we  can 
show  when  estimating  the  radiated  power  that  their  effect  is,  in  general, 
negllg  ible. 

The  reduced  a...  ■  "^al  equations  is  as  follows 
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where  only  the  first  term  of  equation  (3.16)  i b  now  Included,  the  other 
being  negllRihle  as  shown  by  equation  (5.6);  with  a  similar  equation 


for  Y-type  edge  modes: 
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where  SC  c  and  TC  c  are  defined  by  equations  (3,6)  and  (3.19) 
mnqr  mnqr 

respectively:  and 
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The  dominant  terms  in  these  .  educed  equations  have  been  obtained 
from  estimates  based  on  the  numbers  of  resonant  modes.  It  is  easier 
now  to  solve  the  equations  as  they  stand  considering  all  the  modes  in 
the  summations.  The  resonant  modes  will  then  again  oe  picked  out  when 
we  subsequently  average  over  a  narrow  band  of  frequencies. 

Equation  (5.9)  is  almost  identical  to  equation  (4.1),  the  only 
difference  being  that  now  as  n~t  s.1 1  modes  are  corner  modes,  the  summa¬ 
tions  in  equation  (5.10)  are  over  only  a  limited  number  of  modes.  The 
approximate  solution  is  similar  to  equation  (4.11)  with  this  restriction 
applied  to  the  summations. 

The  coupling  terms  in  equation  (5.8^  satisfy  the  relation 
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exactly.  An  exact  solution  is  thus  possible.  This  is 
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A  similar  solution  exists  for  v  ^  .  Finally,  by  substituting  these 

inn 

values  into  equation  (5.10)  we  obtain  the  approximate  result 
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^  -r**ssions  for  the  ;■*  ■'•‘sted  power  spectral  dens..,  can  now  be 
obtained.  We  again  assume  that  equation  (4.15)  applies  to  the  forcing 
field.  The  radiated  power  is  the  sum  cf  the  power  radiated  by  each 
type  of  mode  separately.  We  will,  consider  here  only  the  resistive 
powei.  Thus 


'll  *  +  H  y  +  T\c  -r  Tf r 


These  terms  are  estimated  separately.  We  obtain „  analogous 


to  equation  (4.16)  (where  now  the  real  part  only  is  implied) 
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with  similar  express  ions  for  both  n  (a1)  and  II  (u>).  The  expression 
for  tt^(w)  involves  a  double  summation: 
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This  expression  is  simplified  by  substituting  for  the  coupling  coef¬ 


ficients  and  performing  one  of  the  summations  in  each  of  the  last 

2  2 

two  terms.  The  two  summations  are  over  the  terms  k  and  k  ,  respectively 

n  m 
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Tf  large  numbers  of  acoustically  fast  modes  are  excited  these  summation 

ko£3  ko£l 

have  die  respective  values,  noting  the  parity  requ ireinent „  — and  j  — 

v 

Hie  radiation  excited  by  the  term  <+>  taking  both  edge  and  acoustically 

mn 

fast  modes  into  account  is  then  found  to  be 
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Tne  direct  radiation  at  resonance  is  considerably  higher  than  the 
radiation  from  the  coupled  nonresonant  acoustically  fast  modes.  We 
thus  obtain  as  an  approximation  to  the  radiated  power  spectral  density 
tne  sum  of  the  radiation  from  each  type  ct  mode  independe  -iy. 
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The  summations  are  nominally  over  all  modes  of  each  type.  "nwever, 
'Tc  r-*>!i;i  '.on  of  tiie  modal  equations  to  a  form  we  could  solve  was  made 
o'  the  TaoU  rf  ect  In a ten  of  1»*  between  resonant  modes.  These 

•nodes  are  picked  out  when  equation  (S.ll)  is  averaged  over  a  narrow 
band  of  frequencies.  The  ••  upling  effects,  as  in  the  low  frequency 


rise,  are  contained  in  the  admittance  functions  only.  We  have  gone  to 


some  lengths  to  show  that  most  cf  the  coupling  effects,  In  particular 
those  better*  different  types  of  modes,  ai.e  small.  The  neglect  of  the 


additional  coupling  terms  will  he  most  justified  when  the  structural 
damping  is  small.  No  modal  coupling  now  occurs,  except  between  different 
modes  with  the  same  resonance  frecuencv.  We  do  not  include  this  pos¬ 
sibility  in  the  analysis  The  criteria  to  be.  satisfied  for  there  t; 
be  no  coupling  are: 

i)  no  coupling  of  comer  modes  if 
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ii)  no  coupling  of  edge  modes  if 
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where  and  rre  the  modal  densities  of  those  edge  modes  that  are 
coupled  defined  by 
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When  the  above  criteria  are  satisfied,  averaging  equation  (^.H) 
over  a  narrow  band  of  frequencies  leads  to  the  straightforward  result 
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In  this  expression,  the  edge  mode  damping  and  modal  correlation 

A  3k 

functions  are  treated  as  uniform  for  all  edge  modes,  NX  -  -  --£  k  -r-™—  Aw 

&  r  o  3  uj 

F  A  ko 

is  the  number  of  resonant  edge  modes  of  each  type,  and  N  »  °  is  the 

4  TT 

nuiaber  of  acoustically  fa„t  modes.  This  result  is  the  same  as  pre¬ 
viously  obtained  light  fluid  loading  results  (e.g.  references  2  and  4), 
with  the  addition  of  non re so nan t  acouscically  fast  modes.  Hie  effect 
of  these  additional  modes  is  ,  In  general,  small  as  is  shown  in 
section  7. 

When  the  criteria  (5.121  ^e  satisfied  lde  rah  1  e  modal 

coupling  occurs.  The  frequency  averages  over  the  admittance  functions 
must  now  be  obtained  as  in  section  4  (equation  (4.20),  forward).  We 
obtain  the  result,  analogous  tc-  equation  (4. 30)  in  the  low  rrecuencv  case 
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The  (m,n)  summat tons  and  averages  here  are  over  all  the  resonant 
modes  of  each  type  in  the  ,\u  frequency  nand.  The  (q,r)  average  is  over 
all  the  acoustically  fast  modes. 

The  additional  modal  inertia  terms  arising  from  the  admittance 
functions  again  serve  only  to  modify  slightly  the  modal  resonance 
treq\:ency  and  hence  the  determination  of  which  modes  are  resonant  in 
the  in)  frequency  hand.  Their  effect  Is  again  small;  the  important 
inertia  term  is  the  self  inertia  defined  bv  equation  (3.20). 

Typical  numerical  values  of  equations  ( h . 133  and  (b.14)  are  d  is 
cussed  in  section  l7. 


!■  Pie  intensity  Directivity  cf  the  Radiated  Fie  ■  d 

We  will  discuss  brief Iv  the  directivity  pattern  of  the  inter 
sit”  of  the  radiated  field.  At  suff icientlv  large  distances  free 
plate  the  far  f  f  *  1  d  approximation  co  the  radiated  pressure  ‘ feld 
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the  only  component  being  along  the  radius  vector  from  the  origin. 

We  will  assume  that  equation  (4.15)  again  applies  to  th  •  external 

F 

field  and  neglect  the  effect  of  $  .  We  then  obtain  a  simple  anproxi- 

mn 

mate  form  lor  the  intensity  of  the  acoustic  field: 
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which,  on  averaging  over  a  narrow  frequency  hand  reduces  to 
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The  summation  in  equation  (6.1)  is  over  all  modes  with  resonance 
irequer  'es  in  the  hi)  frequency  band.  The  expression  is  valid  at  all 
(sub-acoustic  critical)  frequencies. 

2  2 

The  directivity  is  determined  by  the  terms  |l  |  and  1 |  • 

At  low  frequencies  the  radiation  obviously  uniform  in  direction. 

At  sufficiently  high  frequencies,  the  dominant  radiation  is  from 

r  example,  an  X-type  edge  mode  (k  <  k  ,  k  >  k  ) 

mono 


edge  modes.  Now, 


has  a  sharp  maximum  intensity  of  radiation  in  the  direction 
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We  can  approximate  j  I„  j  hare  by  2/k  .  For  anv  given  mode,  the 

3n'  mn  '  ° 

direction  of  maximum  radiation  at  its  resonance  frequency  is  nearer  to 

the  normal  to  the  panel  under  fluid  loading  conditions  because  the 

x, 

corresponding  k  is  decreased.  The  change  in  —  -  is  in  the  ratio  of 
O  R 

the  change  in  resonance  frequencies,  namely 


But,  at  any  given  frequency  to,  k^,  and  hence  the  "width"  of  the  edge 
mode  regions  in  k-space ,  is  fixed.  Thus,  for  example,  the  X-mode 
radiation  is  always  due  to  modes  with  the  same  values  of  k  ,  although 

HI 

the  values  of  k^  corresponding  to  the  resonant  modes  at  this  frequency 
will  change.  The  directivity  pattern  of  radiated  intensity  at  high 
frequencies  is  thus  essentially  independent  o  luid  loading,  although 
the  overall  magnitude  of  the  intensity  will  depend  on  changes  in  the 
wavenumbers  associated  with  a  resonance  frequency,  and  on  the  relative, 
numbers  of  resonant  edge  modes  in  a  bandwidth. 
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We  note  finally  that  the  radiated  power  is  the  total  intensity 
integrated  over  a11  directions.  To  demonstrate  this  we  require  from 
equation  (6.1)  the  in.  gral 
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where  Pi  is  the  surfac*  of  solid  angle  2rr  at  radius  K.  For  the 
X-edge  modes  we  make  the  approximation 
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Using  an  approximation  similar  to  that  used  in  section  3, 


namely 
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we  can  obtain: 
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with  a  similar  expression  for  i?.  ^  is  easily  calculated  by  approxi- 
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mating  1 1.  j  and  |l,  j  as  2/k^  and  2/k  ,  respectively,  giving 
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Using  these  values  together  with  equation  (6.1)  we  finally 


obtain 
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This  expression  is  equivalent  to  equation  (5.13)  if  <i>  is  i-ero 
and  when  suitable  expressions  are  substituted  for  the  modal  radiation 
coefficients;  and  to  equation  (A. 29)  when  only  corner  modes  are  excited. 


7 .  Evaluation  o f  the  Effect  of  Fluid  Loading 

In  this  section  some  typical  values  of  the  power  spectral  densities 
of  radiated  sound  obtained  in  sections  A  and  5  are  estimated  when  the 
plate  is  excited  by  a  turbulent  boundary  layer.  Corcos ' ^  modei  is 
used  for  the  cross-spectral  densitv  of  the  pressure  field  generated  by 
the  turbulent  boundary  layer.  This  Is 
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where  we  have  further  assumed  that  the  lcn;gitudinal  end  lateral  ampli¬ 
tudes  decay  exponentially.  In  this  expression  4>(w)  is  the  spectral 
density  of  the  mean  square  boundary  layer  pressure,  =  (r^,r.,)  is  the 
spatial  separation,  is  the  convection  speed  of  the  pressure  field 
and  anda^  are  non-dimens ionar  constants.  In  reference  4  it  is 
shown  that  the  modes  are  uncoupled  (i.e.  condition  -.,15  holds)  if 
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it  is  assumed  that  these  inequalities  are  satisfied  for  the  range 
of  frequencies  we  will  consider.  It  is  further  assumed  that  Uc  is  so 
small  that  no  hydrodynamic  coincidence  effects  occur.  This  last 
assumption,  although  reasonable,  is  maue  solely  to  make  the  deter¬ 
mination  of  d>  easy  and  does  not  affect  the  discussion  o  the  fluid 
mn 

loading  effects.  Under  these  assumptions  we  obtain  the  simple  relation 
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i>  Is  thus  constant  for  all  inodes  being  a  function  only  of  fre- 
rran 

quency . 

The  structural  damping  of  the  plate  is  an  important  factor  in 
determining  the  velocity  response  and  the  associated  radiation.  Un¬ 
fortunately  no  simple  description,  if  any,  of  the  structural  damping  is 
available.  Most  of  the  energy diasipa red  in  thin  plates  is  into  the  end 
supports  and  thus  depends  greatly  on  the  means  of  supr^rt  and  the  struc¬ 
tures  to  which  the  plate  is  attached.  The  ef 'eet  of  the  structural 
damping  is  further  modified  under  dense  fluid  loaning  conditions  by 
its  magnitude  reljtive  to  the  radiation  damping.  The  addition  of  damping 
treatment  to  a  plate  v.  an  lead  to  different  effects  under  light  and  dense 
fluid  loading  conditions.  To  demonstrate  the  interaction  of  both 
structural  and  fluid  damping,  and  the  magnitudes  of  the  fluid  loading 
effects,  some  extremely  simplyfying  assumptions  will  be  made  about  the 
structural  damping.  Ke  will  assume  that  the  structural  loss  factor  is 
inversely  proportional  to  frequency  and  is  independent  of  wavenumber. 

This  last  assumption  ip  equivalent  to  assigning  the  total  structural 
loss  factor  measured  in  narrow  frequency  bands  to  each  mode  in  the  band. 
For  thin  plates  the  total  structural  loss  factor  is  found  in  genor. 


to  be  inversely  proportional  :c  frequency,  but  is  usuallv  measured  on.y 
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over  fairly  sir  'll  ranges  of  frequency.  We  will  partly  compensate  for 
this  lack  of  knowledge  of  values  of  the  loss  factor  by  evaluating  the 
power  for  several  widely  different  values.  We  thus  write  for  all 
>  modes 


I.,  =  1 


A 


where  ji  is  constant,  and  evaluate  the  radiated  power  using  various 
values  of  ^  . 

Under  these  simplyfying  assumptions,  the  evaluation  of  the  radiated 
power  spectral  density  is  particularly  straightforward.  For  example, 
at  low  frequencies,  by  averaging  over  all  modes  in  a  band,  we  obtain 
from  equation  (4.29) 
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and  from  equation  (4.30): 
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The  average  values  of  the  radiation  coefficients  obtained  ix\  reference 


4  are  used.  These  differ  bv  a  factor  2  fro*'  those  published  earlier 


(2) 

by  Maidanik  .  No  account  will  be  taken  of  the  second  order  inertia 
terms  when  alculatlng  the  modal  resonance  frequences.  In  the  calculations 
performed  here ,  the  wavenumber  corresponding  to  resonant  modes  at  any 
frequency  (including  the  modal  9elf- inertia  term),  that  is,  the  solu¬ 
tion  of  the  equation 
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was  obtained  graphically. 

Typical  values  of  these  expressions  are  compared  in  Figure  3  with 
the  spectral  density  obtained  by  neglecting  fluid  loading,  that  is, 
with  the  expression 


where  here  the  modal  expressions  are  evaluated  at  the  in  vacuo  r  onant 
frequencies  live  expressions  obtained  in  section  5  can  be  simplified 
In  a  similar  manner. 

Figures  3  and  6  show  some  typical  values  of  the  power  radiated  by 
a  2'  x  2'  x  1/10"  steel  plate  water  loaded  on  one  side.  A  very  thin 
piate  is  chosen  to  demonstrate  more  markedly  the  effects  of  fluid  load¬ 
ing.  For  thicker  plates  the  modal  density  is  too  low  for  the  averag¬ 
ing  technique  over  resonant  modes  in  narrow  frequency  bands  to  In 
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applicable  at  low  frequencies  (assuming  1/3  octave  band3  are  used)  , 
although  the  radiation  may  be  readily  computed  from  equation  (A. 19). 
Three  different  values  of  6  are  used,  namely,  0  *  tt  ,  20 tt  and  200  tt. 

These  values  correspond  to  quality  factors  at  1000  Hertz  of  2000,  100 
and  10,  respectively.  The  range  thus  covers  both  very  lightly  and  very 
heavily  damped  systems.  We  would  expect  the  value  <J>  *  100  at  1000 
Hertz,  t"  be  of  meet  practical  interest. 

The  curves  shown  in  Figures  5  and  6  are  not  restricted  to  the  low 
and  high  frequency  regimes  considered  in  sections  4  and  5  (in  this 
example,  these  regimes  correspond  to  f  <  625  Hz  and  to  f  >>  1250  Hz). 

Rather  it  is  assumed  that  che  corner  modes  are  always  essentially  un¬ 
coupled  to  other  types  of  modes  and  that  the  corner  mode  solution  is 
thus  applicable  at  all  frequencies.  To  this  solution  must  be  added 

the  edge  mode  radiation  when  k  £„ ,  k  (,  >  3r.  (In  reference  s  it 

oi  o  J 

is  noted  that  the  first  few  edge  modes  excited  as  one  considers 

increasing  frequencies  have  radiation  coefficients  more  characteristic 

of  corner  modes.  ill's  fact  has  been  used  to  estimate  the  un-loaded 

radiation  and  accounts  for  the  discontinuities  in  the  spectra  at 

k  v, ,  k  C.  -  3r.  The  edge  mode  radiation  under  dense  fluid  loading 

o  1  o  ! 

is  shown  dashed  in  Figure  5  for  frequencies  v  <  k  ?  <  3".)  f.xcept 

for  the  very  heavily  damped  case,  the  edge  mode  radiation  is  less  than 

the  corner  mode  radiation  until  manv  edge  modes  are  excited.  This 

fact,  together  with  the  fact  that  the  contribution  from  the  acoustically 
F 

fast  mooes  due  to  0  is  negligible  (in  this  example  the  a.f.  radiation 
mn 
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is  less  than  5 %  of  the  total  radiation  at  all  the  frequencies  considered) 
lends  added  justification  for  the  neglect  of  additional  coupling  effects 
in  the  corner  mode  solution  in  the  middle  range  of  frequencies. 

The  effect  of  an  increase  in  the  structural  damping  on  an  un-ioaded 
plate  is  a  uniform  reduction  in  the  radiated  spectrum.  The  effect  on 
a  fluid-loaded  plate  is  less  straightforward.  We  consider  first  the 
comer  mode  radiation  in  Figure  6.  For  $  *  it,  a  very  lightly  damped 
plate,  the  additional  damping  due  to  fluid  loading  causes  a  marked 
decrease  in  the  radiated  spectrum.  For  greater  damping,  P  *  20r, 
the  additional  effect  of  fluid  loading  is  not  as  marked.  Soth  these 
cases  have  been  evaluated  using  equation  (4.29)  (and  its  analogy  at 
higher  frequencies).  However,  further  increase  of  the  structural  damp¬ 
ing,  P  ”  200 r.  causes  modal  coupling  as  the  resonance  peaks  c'-erlap. 
Equation  (4.30)  is  now  applicable.  Figure  6  shows  that  rather  than  a 
further  decrease  in  the  fluid  leading  effe  t  as  :•  is  increased  from  20v 
to  200 r,  the  modal  coupling  results  in  a  greatly  increased  radiation 
damping  and  lienee  shows  •»  marked  fluid  loading  effect.  Further  increase 
of  ,  of  course,  would  again  result  in  a  decreasing  change  due  to  fluid, 
loading.  This  behavior  is  characteristic  of  corner  mode  radiation  at 
all  frequencies. 

The  radiation  damping  of  edge  modes  is  verv  high.  For  a  2  ’  x  .'  * 
stool  pinto  in  water  the  uncoupled  radiation  damning  of  an  edge  mode 
corresponds  to  a  quality  factor  of  20  at  ?'W0  b'r.  (Th  is  value  is  iepen  • 
dent  on  ‘be  thickness  <  !  the  plate  or.lv  through,  the  chance  of  resonance 
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frequency  due  to  mass  loading.)  Unless  the  structural  damping  is  verv 

high,  the  factor  0  c  o  (urn  n  +  0  c  0  )  '  is  essentially  unitv 
o  o  ittr  p  mn  o  o  mn 

for  all  edge  modes.  The  increased  radiating  efficiency  over  the  corner 


modes  is  thus  offset  by  the  increased  damping  and,  as  seen  in  Figure  5, 
verv  manv  edge  modes  must  be  excited  before  the  edge  mode  radiation 


dominates  the  cornier  mode  radiation.  This  is  true  for 


and  f  =  20’ 


At  verv  high  levels  of  damping  however,  r  -  200  r,  the  increased 

corner  mode  radiation  damping  Hue  to  modal  c.-apl  fnc  decreases  th cerr»r 
mode  rad  iat ion  to  such  an  extent  that  edge  mode  rad iat ion  again  alvavs 
predominates.  Th  large  discont  inuff  in  the  curves  for  -  -  2>V1 ~ 

presumable  '.;e  to  our  lnahilit"  to  anal  vise  compietolv  the  modal  coup¬ 
ling  effects  in  the  middle  .Hina-  of  f  cecuenc  ie-< .  As  las  '-eott  nrev  ious  1  v 
remarVed  the  analysis  given  in  section  •->  in  5  "  is  most  accurate  for 
lev  valves  of  tl’  structural  damning. 


done  1  us  ion 


he  modal  avnrcav1-  used  her*'  is  perhaps  not  part  !«.  t  isriv  suited 
;  t  ••  s . ;  v  of  the  .>f  foots  of  fluid  loading  as  it  leads  to  an  Inf  in  ite 
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At  low  frequencies,  k  £  <  r'^ '  t^‘e  aPProxlmations  t0  the 

coupling  coefficients  of  the  modal  equations  are  good  and  the  resulting 
equations  can  be  solved  fairly  accurately.  The  radiated  power  spectral 
density  can  be  obtained  either  by  averaging  over  resonant  modes  in 
narrow  frequency  bands  when  the  modal  ity  is  high  (equation  (4.20) 
and  (4x30))  or  computed  directly  from  equation  (4.18). 

At  higher  frequencies,  the  complexities  of  the  modal  interactions 
are  such  that  considerable  simplification  of  the  modal  equations  must 
be  made  before  even  an  approximate  solution  is  found.  Hie  approx ir-rt ions 
made  in  this  case  are  more  justifiable  in  cases  where  the  structural 
damping  is  so  low  that  in  any  case  a  negligibl  amount  of  modal  inter¬ 
action  occurs. 

The  main  effects  of  fluid  loading  have  been  discussed  in  section  7. 
We  have  noted  that  the  inertia  coupling  terms  pla)  a  somewhat  minor  role. 
This  is  ot  really  surprising:  there  can  be  no  exchange  of  energy  via 
inertial  coupling,  we  have  treated  systems  such  that  the  response  of 
the  system  at  any  frequency  is  described  by  the  response  of  those  modes 
that  are  resonant  at,  or  near,  that  frequency.  The  response  is  thus 
pri  rily  determined  by  the  amplitudes  of  the  modal  resonances.  These 
amplitudes  are  changed  from  the  in  vacuo  case  solely  by  the  additional 
energy  loss  by  acou;."ic.  radiation  to  infinity:  we  call  this  additional 
energy  loss  the  radiatir”  damping  of  the  plate.  The  slight  changes  in 
resonance  frequencies  caused  by  the  inertia  coupling  are  overshadowed 
by  the  large  changes  in  frequency  caused  by  the  modal  self  inertia. 
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The  magnitude  of  the  radirtion  damping  is  partly  deter., ined  by  the 
amount  of  modal  interaction  that  occurs.  When  the  structural  damping 
is  small  the  resonance  peaks  are  very  sharp;  no  modal  interactions  occur 
as  the  energy  of  the  system  is  contained  in  very  narrow,  separate  bands 
of  frequencies.  Thus,  the  acoustic  fJeld  at  any  frequency  is  generated 
solely  by  the  mode  that  is  resonant  at  that  frequency.  The  radiation 
damping  of  each  mode  can  then  only  be  due  to  the  acoustic  field  generated 
bv  that  mode  alone.  Except  for  very7  thin  plates,  this  is  the  most 
typical  situation  met  with  in  practice.  When  the  structural  damping 
is  large,  the  widths  of  the  resonance  peaks  are  increased,  i.e. ,  each 
mode  is  considerably  excited  over  a  wider  band  of  frequencies.  The  acous¬ 
tic  field  at  any  frequency  is  now  due  not  only  to  the  mode  resonant  at 
that  frequency  but  also  to  other  modes  with  resonance  frequencies  near  to 
that  frequency.  The  radiation  damping  is  correspondingly  hi&ner.  The 
net  result,  of  course,  is  a  decrease  in  the  total  radiated  field  because 
of  the  decreased  amplitude  of  the  plate  response. 

The  two  important  features  of  fluid  loading  are  the  fluid 
inertia  effect  and  the  radiation  damping.  Both  effects  lead  to  a 
decrease  in  the  radiated  acoustic  field.  The  radiation  damping  causes 
a  decrease  in  the  velocity  response  amplitude  of  the  plate.  The  fluid 
inertia  causes  o  decrease  in  the  modal  resonance  frequencies.  At  any 
f'teq,.  ncy  the  resonar.  „<des  correspond  to  higher  wavenumbers,  and  thus, 
have  lower  radiation  efficiencies. 
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APPENDIX 

It  has  been  assumed  throughout  the  analysis  that  condition  4.15 
applies  to  the  external  force,  that  is,  there  is  no  modal  coupling 
induced  by  ^he  applied  force.  However,  at  low  frequencies  a  simple 
expression  for  the  radiated  power  similar  to  equation  (4.18)  can  still 
be  obtained  even  if  this  condition  does  not  hold.  The  many  cross-coup¬ 
ling  terms  representing  power  flow  between  the  modes  again  cancel  and 
we  obtain  the  expression 
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where  we  have  made  use  of  the  approximate  result  4.14. 

The  spectrum  of  the  radiated  power  is  represented  by  the  veal 
part  of  equation  (A.l).  This  expression  can  be  averaged  over  a  narrow 
band  of  frequencies  in  cases  where  the  modal  density  is  high.  We  require 
an  estimate  of  the  integral 
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Using  the  notation  of  section  4,  and  assuming  light  structural 


damping,  this  Integral  has  the  form 
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This  result  is  a  factor 
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2 

times  the  corresponding  integral  over  |y  |  .  We  assume  that  n  /e, 

mn  mn  X 

is  always  very  small,  that  is,  the  widths  of  the  modal  resonance  peaks 
arc-  less  than  the  width  of  the  frequency  band  over  which  we  Integrate. 
The  corner  mode  results  given  in  section  4  will  thus  be  approximately 
correct  for  low  values  of  the  structural  damping  even  when  condition 
4.15  does  not  hold.  For  excitation  by  a  turbulent  boundary  layer  the 
conditions  for  equation  (4.15)  to  hold  are  frequency  dependent.  It  is 
at  low  frequencies  that  the  additional  approximation  give*  here  may 
be  applicable. 
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